Ferromagnetic coupling of mononuclear Fe centers in a self-assembled metal-organic 

network on Au(lll) 



T. R. Umbach\ M. Bernien^ C. F. Hermanns^ A. Kriiger^ V. Sessi^, 1. 

Fernandez-Torrente\ P. Stoll^'^, J. 1. Pascual^'^'"^, K. J. Franke^'^'^ W. Kuch^'^ 
^ Freie Universitdt Berlin, Fachbereich Physik, Arnimallee I4, 14-195 Berlin, Germany 
^ European Synchrotron Radiation Facility, PB 220, 38043 Grenoble, France 
^ Genter for Supramolecular Interactions, Freie Universitdt Berlin, Arnimallee I4, 14195 Berlin, Germany 
^ GIG nanoGUNE, 20018 Donostia-San Sebastian, and ikerbasque, 
Basque Foundation for Science, 4^011 Bilbao, Spain 
^ Institut fiir Festkorperphysik, Technische Universitdt Berlin, Hardenbergstrafie 36, 10623 Berlin, Germany 

(Dated: December 17, 2012) 

The magnetic state and magnetic coupling of individual atoms in nanoscale structures relies on a 
delicate balance between different interactions with the atomic-scale surrounding. Using scanning 
tunneling microscopy, we resolve the self-assembled formation of highly ordered bilayer structures of 
Fe atoms and organic linker molecules (T4PT) when deposited on a Au(lll) surface. The Fe atoms 
are encaged in a three-dimensional coordination motif by three T4PT molecules in the surface plane 
and an additional T4PT unit on top. Within this crystal field, the Fe atoms retain a magnetic ground 
state with easy-axis anisotropy, as evidenced by X-ray absorption spectroscopy and X-ray magnetic 
circular dichroism. The magnetization curves reveal the existence of ferromagnetic coupling between 
the Fe centers. 



PACS numbers: 75.70.-i, 78.70.Dm, 68.37.Ef, 75.30.Et 

As the size of a magnetic material decreases, the sta- 
bility of the collective magnetic groundstate becomes 
weaker, leading to a paramagnetic behaviour at atomic 
scale dimensions [1]. The identification and study of 
mechanisms for stabilizing magnetically oriented ground 
states aims to provide routes for the design and construc- 
tion of nanomaterials (clusters and thin films) that retain 
some magnetic order at temperatures as high as possible. 
At the limit of individual atoms deposited on a metal sur- 
face, magnetic anisotropy and interatomic coupling are 
the key ingredients to tune the magnetism of a nanostruc- 
ture 0. The anisotropy of the atomic environment may 
lower the magnetic moment but, in turn, defines easy 
directions for magnetization. Exchange interactions cou- 
pling the magnetic moment of neighbor atoms drive the 
formation of magnetically ordered states, but this effect 
quickly vanishes with increasing interatomic distance. 
Only the indirect exchange coupling through valence 
electrons of the metal support (the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction [3U7]) may be op- 
erative. It has been shown that RKKY interactions lead 
to collective ferromagnetic (FM) or anti-ferromagnetic 
(AFM) ground states, but the exchange energy is small 
(J -50 /ieV) EH [8]. 

An interesting approach to enhance the magnetic 
anisotropy of atoms in a thin film is to augment the ligand 
field splitting of the d-levels by coordination to specifi- 
cally selected organic ligands [2 . The formation of a co- 
ordination bond to molecular endgroups induces a charge 
redistribution at the metal sites and cast the ligand field 
by adopting high-symmetry metal-organic structures and 
periodic networks O |9l HO] . In bulk metal-organic mate- 
rials, the organic linkers connecting magnetic metal cen- 



ters also mediate the magnetic coupling between them 
by means of electronic states extended along the organic 
back-bone [T2j [13] . The presence of such superexchange 
mechanism for magnetic coupling in metal-organic thin 
films could optimize the stabilization of collective mag- 
netic ground states, but its detection remains still a chal- 
lenge. 

Here, we show the existence of ferromagnetic coupling 
between Fe ions within a self-assembled metal-organic 
film grown on a Au(lll) surface. Single Fe atoms lie in 
three-dimensional coordination sites and are connected 
to neighboring Fe atoms by 2,4,6-tris(4-pyridyl)-l,3,5- 
triazine (T4PT) molecules (see Fig. ^p)- Analysis of 
X-ray absorption spectroscopy (XAS) and X-ray mag- 
netic circular dichroism (XMCD) measurements reveals 
a high spin S=2 ground state with an out-of-plane mag- 
netic anisotropy. Magnetization curves of the ultra-thin 
layer exhibit deviations from the case of non-interacting 
paramagnetic moments, which are interpreted as a finite 
ferromagnetic interaction between the mononuclear cen- 
ters on the basis of a mean-field model. 

Fe-T4PT networks were grown on an atomically clean 
Au(lll) single crystal. The gold surface was first pre- 
covered at room temperature by a submonolayer amount 
of T4PT molecules, sublimated from a Knudsen cell, and 
by Fe atoms. Subsequently, the substrate was annealed 
to 400 K to activate the formation of coordination bonds. 
The structure of the resulting surface was characterized 
in situ, in our low temperature (4.8 K) scanning tunnel- 
ing microscope (STM). XAS and XMCD measurements 
in total electron yield were performed at the beamline 
ID08 at the European Synchrotron Radiation Facility 
(ESRF). The beamsize on the sample was about 1.4 mm 
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FIG. 1: (a) Large scale STM topography image of the coor- 
dinated Fe-T4PT islands on the Au(lll) surface (/= 253 pA, 
Vs — —1.5 V). The Fe-T4PT islands are homogeneous over 
the whole surface. There is no indication of Fe clusters. The 
STM image was obtained with a room temperature STM 
at the ID08 beamline at the European Synchrotron Radi- 
ation Facility (ESRF) directly before X-ray investigations, 
(b) Molecular structure of 2,4,6-Tris(4-pyridyl)-l,3,5-triazine 
(T4PT). (c) An STM zoom evidences the metal coordination 
of the T4PT molecules (/= 230 pA, Vs = 0.06 V). STM image 
processing by WSxM ^4]. 



(FWHM) in the horizontal, and 50 /im in the vertical 
direction. The sample was kept at a temperature of 8 K 
in a variable external magnetic field up to 5 T and ex- 
posed to X-rays under different incidence angles. Here, 
STM measurements (Fig. [T]) were carried out in-situ at 
several macroscopically separated areas, indicating that 
the sample was covered by a homogeneous metal-organic 
network, and that there are no uncoordinated Fe clusters. 

Large-scale STM images show the tendency of Fe 
atoms and T4PT molecules to coordinate in networks 
with high a degree of order (Fig. [l^). The presence 
of the Au(lll) herringbone reconstruction underneath 
the metal-organic layer indicates a weak interaction with 
the substrate. STM images recorded at low bias voltage 
(Ks < 0.5 V) resolve the skeleton of the T4PT molecules 
(Fig. [l]^). Each pyridyl group participates in a threefold 
node, a configuration that can only be stable if the natu- 
ral repulsion between electrophilic endgroups is overcome 
by insertion of a bonding center prone to coordinate with 
the N lone-pair electrons |2l|9j. Hence, we anticipate the 
inclusion of Fe atoms at these sites [15] . 

The structure of the Fe-T4PT network shows a pecu- 
liar dependence on the bias applied during its inspection 
(Figj2|i and b). Starting from the low bias structure in 
Fig.|2p!, an increase of the sample bias above 0.6 V results 
in the observation of a new set of T4PT molecules cen- 
tered on the Fe sites and rotated by = 29° (Fig. ^p). 
These new features can not be associated to tunneling 
through different orbitals of the same molecules due to 
their different orientation. Furthermore, we observe do- 
mains with mirror symmetry around dislocation lines at 
high bias voltage that are absent when inspected at low 
bias voltage (see supplementary information SI). Based 




FIG. 2: (a)-(b) STM images of a Fe-T4PT network measured 
at different sample bias voltages (/= 0.33 nA, (a) Vs = 0.31 V, 
(b) Vs = 2.0 V). (c) Structure of the first-layer of the Fe- 
T4PT network. Fe atoms are three- fold coordinated by three 
T4PT molecules in the first layer, (d) In the second layer, 
T4PT molecules are centered at the Fe sites and are rotated 
by = 29° with respect to the first-layer molecules (indicated 
by the black lines). 



on these facts, we conclude that the film consists of a 
bilayer of T4PT molecules: Fe atoms coordinate to the 
pyridyl groups of three T4PT in the bottom layer (first 
layer in Fig. ^) and to the triazine center of an addi- 
tional T4PT molecule directly on top (second layer in 
Fig. [2]i). Their different aspect with the bias voltage 
can be rationalized by considering that the top layer is 
only weakly coupled to the metallic substrate and poorly 
screened [16]. Its electronic gap becomes wider and the 
molecule more transparent to electrons tunneling close 
to the Fermi level Ep [l7l|18], allowing us to resolve the 
structure of the bottom layer at low bias voltages. 

In the resulting structure the Fe atom is located in a 
three-dimensional coordination cavity, encaged by a to- 
tal of six nitrogen atoms (Fig. ^ and d). The presence 
of the top T4PT centered directly on Fe sites suggests 
that the triazine moiety also participates in the coordi- 
nation with the metal atom. A probable scenario is that 
this metal-ligand bond causes an upward displacement 
of the Fe atom [2l [191 1201 ? which would cause a distor- 
tion of the ligand field around the Fe atom, bringing it 
from a trigonal-planar into a three-dimensional (pyrami- 
dal) structure [211 [22] . This configuration could help to 
prevent the quenching of its magnetic moment due to the 
interaction with the gold surface. 

To investigate the magnetic properties of the Fe atoms 
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FIG. 3: (a) X-ray absorption spectra and (b) XMCD dif- 
ference spectra at the Fe L2,3 edges of Fe-T4PT on Au(lll), 
measured at a temperature of 8 K in an applied magnetic 
field of 5 T for two different angles of incidence, as defined in 
the sketch in the top panel. The magnetic field was always 
parallel to the X-ray beam. The dashed line represents an 
XAS of the L2,3-edge in Fe2 03, where Fe is in a +3 oxidation 
state. The shift of the absorption edge is due to the chemical 
sensitivity of the binding energy of the core-level electrons. 



in the Fe-T4PT networks, we have performed XAS and 
XMCD measurements at T = 8 K in an applied magnetic 
field of 5 T (Fig. [3|. From the line shape and position 
of the XAS Ls peak (Fig. [3|i, maximum at 707.2 eV) 
we conclude that we probe isolated Fe atoms in a +2 
oxidation state [23l f2&-28 (the spectrum clearly devi- 
ates from that of Fe(III), shown with a dashed line in 
Fig. 3a for comparison). A +2 oxidation state is typical 
in bulk mononuclear iron complexes [22 , thus support- 
ing the three-dimensional coordination motif described in 
Fig. [2]i. Multiplet calculations of the XAS and XMCD 
difference spectra indicate an = 2 high-spin state (see 
supplementary information S2). 

The presence of a magnetic moment is further corrob- 
orated by XMCD spectra shown in Fig. 3b. The XMCD 
difference signal of the L3 peak exhibits a characteristic 
dip-peak structure at the low-energy side. While the line 
shapes are similar for different angles of incidence, there 
is a large change in the XMCD peak height: the signal 
at normal incidence is more than a factor of two higher 
than at grazing incidence (20°). The higher XMCD sig- 
nal in normal incidence points towards an easy-axis mag- 
netic anisotropy along the direction normal to the sur- 
face, which is a result of the spin-orbit-induced mixing of 
nearly degenerate Fe 3d levels in the surrounding ligand 
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FIG. 4: Magnetic field dependence of the XMCD at the Fe 
Ls edge of Fe-T4PT on Au(lll) for X-ray incidence normal 
to the surface (circles) and at 20° grazing incidence (squares), 
measured at a temperature of 8 K. The magnetic field was al- 
ways parallel to the X-ray beam. Dashed lines are the best 
fits of a paramagnetic spin Hamiltonian for 5'=2, T=8 K, and 
D = —1.16 meV, solid lines are the result of a simulation us- 
ing mean- field-coupled magnetic moments with Tc = 1.8 K 
(see text). The inset shows the mean squared deviation be- 
tween the experimental field dependence of XMCD data and 
the simulation by the paramagnetic spin Hamiltonian (dashed 
line) and the mean- field model, as described in the main text, 
as a function of the anisotropy, represented by the zero-field 
splitting parameter D. For each value of D, all other fit pa- 
rameters have been optimized. 



field. 

Fig. |4] shows the evolution of the Fe L3 XMCD 
signal with applied external magnetic field for both 
perpendicular and grazing incidence. The curvature of 
these plots allows us to quantify the anisotropy of the 
system. We start first by simulating the magnetization 
curves assuming non-interacting paramagnetic moments 
under the presence of magnetic anisotropy by means of 
a spin Hamiltonian approach: 



n = DS^,-I^B9B-S (1) 
The first term describes the anisotropy energy through 
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the zero-field-splitting parameter D and S'^, the spin 
component normal to the molecular plane (a negative 
D stands for an easy axis along z). The second term ac- 
counts for the Zeeman energy (S is the spin vector, B the 
vector of the external field, jj^b Bohr's magneton, and we 
consider a g factor of 2). Solving Eq. [l]for the two direc- 
tions of measurement and calculating the magnetization 
from the thermal population of the resulting levels al- 
lows us to simulate magnetization plots like in Fig. |4] for 
several conditions. The curvature of the magnetization 
plots is an expression of the anisotropy of the system 
and can be tuned by adjusting the parameter D. The 
larger curvature in normal incidence corroborates that 
this direction is the easy axis for magnetization {D < 0). 
However, the experimental curvature cannot be properly 
reproduced, no matter what value of D is used. The 
dashed line in the inset of Fig. [4] shows the mean squared 
deviation between experiment and simulations evaluat- 
ing Eq. [l] as a function of D. We obtain large standard 
deviations for all values of even when the curvature 
is simulated using unrealistically high D values. We thus 
conclude that the experimental data are not compatible 
with non-interacting anisotropic paramagnetic moments, 
but bear instead evidence for a magnetic coupling be- 
tween neighboring Fe moments [29] . 

To include magnetic interactions between nearest 
neighbors, we resort to a mean-field description to sim- 
ulate the magnetization curves. In this model, the mag- 
netic interaction of a spin with all other spins in the 
neighborhood is replaced by an effective mean magnetic 
field added to the external applied field and proportional 
to the magnetization of the ensemble with proportional- 
ity factor Tc=^ Ji/(3/c^). Here Ji is a Heisenberg-type 

i 

interaction of the spin with a spin at site i, is the 
Boltzmann constant. The sum goes over all sites with 
non- vanishing J^. Tc represents the Curie temperature 
in the limit of vanishing anisotropy. 

The solid lines in Fig. [4] are the result of such a simula- 
tion with D = -0.35 meV and Tc = 1.8 K. These values 
correspond to a global minimum (see solid line in the in- 
set of Fig. 4 and supplementary information S4), repro- 
ducing correctly the experimental magnetization curves. 
While the model may not capture the details of the cou- 
pling mechanism, it clearly evidences the necessity to in- 
clude ferromagnetic coupling of the Fe atoms to explain 
the experimental curves. If we assume a coupling to only 
the six nearest neighbor Fe atoms, the coupling constant 
Ji = ks Tc/2 has to be of the order ^ 80 /ieV to repro- 
duce the experimentally observed magnetization curves 
[31]. Since the measurements are carried out well above 
Tc, the effect of correlated fluctuations, which would lead 
to an underestimation of the derived coupling strength in 
the mean-field description, are small. 

Due to the relatively large distance between the Fe 
centers (r = 1.3 nm) a direct exchange coupling is 



very unlikely to be the origin of the observed ferro- 
magnetic coupling. An indirect coupling mechanism 
via the polarization of the surface and bulk conduc- 
tion electrons has been shown to exist on metal surfaces 
[61 [71 [321133]. The RKKY-interaction exhibits an os- 
cillating ferro- / anti- ferromagnetic coupling with the dis- 
tance r of the magnetic moments. It can be described as 
Jrkky ^^(^2^^rY IS] 5 where /ci? is the Fermi wave 
vector, and d defines the dimensionality of the conduction 
electron system. The Fermi wave vector of the Rashba- 
split surface state amounts to kp^i ^ 1.7 nm [35 . Hence, 
a pairwise interaction of the Fe atoms at a distance of 
1.3 nm via RKKY interactions would result in a weakly 
coupled anti-ferromagnetic ground state. Furthermore, 
we expect that the sixfold ligated Fe atoms are lifted 
from the surface [I9l |20] and thus only interact weakly 
with the substrate electrons. Despite small effects that 
may change the scattering phase shift of the conduction 
electrons, we may tentatively exclude a significant con- 
tribution of an RKKY-mediated coupling. 

The molecular ligands coordinated to two neighbor Fe 
atoms can also play an important role in mediating mag- 
netic interactions through a superexchange mechanism. 
Ferromagnetic and anti-ferromagnetic exchange interac- 
tions in metal-organic frameworks are well known in bulk 
materials [12] for nitrogen-based organic linkers and tran- 
sition metal atoms. The superexchange mechanism is 
enhanced when conjugated molecular states bridge the 
magnetic centers [36], as it is the case for Fe-T4PT net- 
works. The spin polarization qualitatively follows a sim- 
ple alternation rule: an even number of linker atoms in 
the pathway leads to anti-ferromagnetic coupling, while 
an odd number results in ferromagnetic coupling. In the 
case of the Fe-T4PT network, the alternation rule results 
in a ferromagnetic coupling between neighboring Fe cen- 
ters. This is consistent with the positive (ferromagnetic) 
Heisenberg-type interaction Ji deduced from the magne- 
tization curves. 

Our results clearly sustain the presence of ferro- 
magnetic coupling of Fe atoms connected by organic 
molecules in a highly ordered metal-organic network. 
The peculiarity of the observed network is that it con- 
sists of a bilayer structure, in which the Fe is coordi- 
nated by three pyridin and one triazine unit. The cen- 
tral triazine ligand may further be responsible for a lift- 
ing of the Fe atoms from the metallic surface due to the 
trans-effect [191 ED!- This configuration eases the for- 
mation of a +2 oxidation state with a net spin. The 
resulting magnetic moments of the Fe ions are coupled 
ferromagnetically with each other, probably via a su- 
perexchange mechanism mediated by the organic link- 
ers. These results further demonstrate that magnetism 
in two-dimensional metal-organic networks on a metal 
surface is feasible without the use of ferromagnetic sub- 
strates and suggests a route to engineer its magnetic char- 
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acter by enclosing the magnetic centers in self- assembled 
three-dimensional coordination cavities. To increase the 
stability of such magnetic networks for potential applica- 
tions, a larger anisotropy and significantly stronger mag- 
netic coupling strengths are required. Chemical design 
strategies towards this goal may comprise the use of rad- 
ical ligands, known from bulk metal-organic ferromag- 
netic materials [371 EH] • 
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STM BIAS DEPENDENCE OF THE FE-T4PT 
COORDINATION NETWORK 

The STM images of the Fe-T4PT network show a 
strong dependence on the applied bias voltage (Figs, [s] 
(a-c)). At low sample bias voltage (0.1 V), we resolve the 
shape of the T4PT molecules, where each pyridyl group 
points to a threefold node. This configuration is stabi- 
lized by the coordination to an Fe atom. At larger sample 
bias voltage (1.0 V in Fig. [s] (b)) the appearance of the 
network changes drastically. We observe triangles formed 
by three distinct protrusions. The center of these trian- 
gular shapes are located on top of the Fe atoms and their 
rotational orientation does not match any high-symmetry 
of the molecular shapes found in Fig. |5](a). Furthermore, 
we note that the orientation of the triangles is not the 
same throughout the whole image (indicated by the white 
line), whereas this is the case for the molecules seen at 
low bias voltage. Hence, we interpret the features at 1 V 
as a second layer of molecules, which may adopt differ- 
ent orientations with respect to the underlying molecular 
layer. The protrusions represent the pyridyl moieties of 
the T4PT molecules. Increasing the sample bias voltage 
further to 3.5 V leads again to a significant change of 
the contrast (shown in Fig. [5](c)). Bright protrusions at 




FIG. 5: (a-c) STM topography images of the Fe-T4PT net- 
work at different sample bias voltages Vs (0.1 V, 1.0 V, and 
3.5 V). At 0.1 V the first molecular layer is visible. At 1.0 V 
the second molecular layer is imaged. The white dashed 
lines indicate a dislocation line in the second layer, where the 
second-layer T4PT molecules change their orientation. This 
change cannot be observed in the first layer (see panel (a)). 
The inset in panel (c) shows a dl/dV spectrum at an Fe center 
with resonances at 0.5 V, 2.0 V, 2.8 V, and 3.5 V. The reso- 
nance at 3.5 V can be associated with a localized electronic 
state of the Fe atom. 



the Fe sites of the Fe-T4PT network are observed. Their 
identification as Fe atoms is sustained by a pronounced 
resonance in the local differential conductance spectra at 
3.5 V at these sites (inset in Fig. [s] (c)), probably aris- 
ing from the coordinated d-states [1, 2 . The appearance 
and disappearance of intensity in the STM images can 
be understood by tunneling through localized molecular 
electronic states. In particular, the second layer is weakly 
coupled to the substrate, thus exhibiting a poor conduc- 
tivity at low bias voltages. The peculiar bias dependence 
is a proof of the two-layer structure of the Fe-T4PT net- 
work: each Fe atom is surrounded by three pyridyl groups 
of three first-layer T4PT molecules and one triazine ring 
of a second-layer T4PT molecule. 



MULTIPLET CALCULATIONS 

The x-ray absorption spectrum and XMCD difference 
spectrum for normal incidence have been simulated by 
multiplet calculations using Cowan's code [3j and the 
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FIG. 6: (a) Helicity- averaged x-ray absorption spectra and 
(b) XMCD difference curves at the Fe L2,3 edges of Fe-T4PT 
on Au(lll). Thick Unes: Experimental spectra, measured at 
a temperature of 8 K in an appUed magnetic field of 5 T in 
normal incidence of circularly polarized radiation. Thin lines: 
Result of a multiplet simulation, see text. 



CTM4XAS program version 5.5 [4 . A good agreement 
for both, the spectral shape of the helicity-averaged ab- 
sorption spectrum and the XMCD difference curve, is 
found for the following crystal field parameters: lODq = 
940 meV, Dr 40 meV, and Da -200 meV. For 
these parameters, which describe a trigonal (Dsd) crys- 
tal field, the Fe ion is in an 6* = 2 high-spin state. Fig. |6] 
shows the resulting simulations (thin lines) together with 
the experimental spectra (thick lines, as also shown in 
Fig. 3 of the main text) for the helicity-averaged x-ray 
absorption spectrum (a) and the corresponding XMCD 
difference curve (b). 



FIT OF THE MAGNETIZATION CURVES 

The field- and angle-dependent XMCD signals pre- 
sented in Fig. 4 of the main text are not directly pro- 
portional to the magnetization M of the Fe ions that is 
modeled by the spin-Hamiltonian. If the magnetic field B 




FIG. 7: Relative alignment of the magnetic field B, the 
magnetization M, and the k vector of the x-rays. 



mean squared deviation (10'^) 




FIG. 8: Two-dimensional color plot of the mean squared de- 
viation of the simulation calculated by the mean-field model 
described in the main text from the experimental field- 
dependence of XMCD data as a function of the magnetic 
anisotropy, represented by the zero-field splitting parameter 
D, and the magnetic coupling between the Fe ions, expressed 
by Tc- The deviation is calculated from the simultaneous 
comparison to the experimental XMCD values for normal and 
20° grazing incidence. The white regions outside the curved 
valley of minimal deviation exhibit deviation values larger 
than 10~^. 



is not applied along the easy axis of the ions, their mag- 
netic moments are not aligned parallel with the magnetic 
field, and the XMCD signal detects only the projection 
of the magnetization onto the k vector of the x-rays. The 
relative alignment of the magnetic field, the magnetiza- 
tion, and the k vector as well as the definition of the cor- 
responding angles are given in Fig. [?[ The XMCD signal 
depends further on the magnetic dipole operator and its 
angle dependence. This term accounts for the anisotropy 
of the spin-density distribution around the Fe ion and can 
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FIG. 9: Magnetic field dependence of tfie XMCD at tfie Fe 
Ls edge of Fe-T4PT on Au(lll) for x-ray incidence normal 
to the surface (circles) and at 20° grazing incidence (squares), 
measured at a temperature of 8 K. The solid black line is the 
best fit of mean- field-coupled magnetic moments with 
the red dotted line for moments with S—1. 



be large for isolated ions. Its temperature dependence is 
the same as the one of the magnetization, and its angle 
dependence is given by (1 — 3cos(26>m)) if the spin-orbit 
coupling is small compared to the ligand-field splitting 
[5]. The intensity of the XMCD signal is calculated from 
the spin-Hamiltonian model according to: 

IxMCD = Acos{Om - Ok) |M| (1 - C(l - 3cos(2^m))) 

where A is an overall scaling factor to the experimental 
signal and C is a factor determining the relative contri- 
bution of the magnetic dipole operator. Besides Tc and 
also A and C were treated as fit parameters for fitting 
the model to the experimental data. 

ERROR MARGINS 

The mean squared deviation between the experimental 
data points of the magnetic field dependence of the 



XMCD difference spectra at the Fe L3 edge and the fit 
regarding to the mean-field model of coupled anisotropic 
magnetic moments is shown as a two-dimensional plot in 
Fig. [8] as a function of the zero-field splitting parameter 
D and the coupling strength, expressed by Tc. In 
addition, C and A have been fitted for each set of D 
and Tc. From Fig. [8] one can see that, although there is 
some correlation between D and Tc, a clear minimum 
is found at the parameters mentioned in the main text, 
i.e., D = -0.35 meV and Tc = 1.8 K. 



SPIN STATE 



Fig. |9] shows a comparison of the best fit of the mean- 
field model described in the main text for = 1 (red 
dashed lines) to the field-dependence of the experimental 
XMCD data (solid symbols) together with the best fit 
for S=2 (black lines). (The latter two are also shown in 
Fig. 4 of the main text.) Similarly good fits are obtained 
for = 1 and S = 2. For 5* = 1 a larger anisotropy and a 
larger coupling between Fe ions result from the fit, with 
values ofD = -1.27 meV and Tc = 3.7 K. 
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